Journal of Genetics and Genomics

While life-history strategies are typically fixed within species, evolutionary transitions between perenniality and annuality can occur. In clonal seagrasses, annual and perennial plants often coexist in the same population, providing a unique model for studying the genetic basis of this transition. Two seagrass Zostera marina populations in South Korea display a striking dichotomy: shallow-water sub-populations follow a typical perennial strategy, whereas their deep-water counterparts are annual. Here we show that this shift from perenniality to annuality, potentially caused by the SAPK7 gene, is genetically coupled with the CAO gene, which is under strong positive selection for low-light adaptation. The up-regulation of the SAPK7 gene triggers early flowering in seedlings, before the formation of any lateral shoots via asexual reproduction. In this special case where the genet contains only one ramet, the post-reproductive death of the ramet is equivalent to the death of the whole genet, which explains the annual phenotype. Our findings reveal a mechanistic example where annuality overcomes perenniality by hitchhiking on a positively selected gene. Given that the ancestral state of plants is perennial, this coupling of annuality with beneficial alleles may represent one of the pathways for the repeated evolution of annual life histories across flowering plants.

Histone proteins, which reside in the nuclei of eukaryotic cells, are involved in diverse cellular processes. The core histone H4 serves as a structural component of the nucleosome. Patients carrying mutations in H4Clustered histone (H4C) genes exhibit a broad spectrum of developmental abnormalities, including short stature, microcephaly, intellectual disability, growth retardation, and digital anomalies. However, the impact of H4 mutations on mammalian embryogenesis remains largely unclear. Here, we demonstrate that histone H4C genes play crucial roles in skeletal development and cortical neurogenesis. We found that mRNAs of the histone H4C gene family are specifically expressed in proliferating progenitor cells in the developing mouse neocortex and in human induced pluripotent stem cell-derived cortical organoids. CRISPR-mediated disruption of H4C3 in mice caused severe defects in skeletal formation and neocortical neurogenesis. Furthermore, overexpression of a mutant form of H4C3 resulted in altered expression of genes associated with cellular migration and motility. Together, these findings suggest that histone H4 plays a critical role in regulating the balance between proliferation and differentiation during mammalian embryonic development, thereby explaining the broad spectrum of patient phenotypes.

The tunicate Ciona robusta offers a tractable non-vertebrate chordate model for probing gene function via tissue-specific, CRISPR/Cas9-mediated mutagenesis in F0. Building on Arcadia Sciences Zoogle platform, which identifies and ranks orthologs of human genes from various non-traditional model organisms, we carried out a pilot project to probe the developmental roles of three notochord- and endoderm-expressed candidate orthologs of human disease genes (Fcho, Pgm3, and Nckap1) alongside a fourth gene (Plastin) implicated in papilla cell elongation. This preprint compiles and updates a series of research project milestones previously posted episodically on Zenodo. Here we summarize the full results and our conclusion about this pilot project. Using CRISPR/Cas9, we found that tissue-specific knockout of Pgm3 and, to a lesser extent, Fcho caused significant defects in larval tail elongation. Separately, CRISPR knockout of Plastin, an actin-bundling gene expressed throughout the sensory-adhesive papillae of the larva, caused a subtle reduction in papilla cell elongation when combined as a duoble knockout with another actin-bundling protein-encoding gene, Villin. These results identify Pgm3 as the most promising candidate for further development as a Ciona-based model of human disease and demonstrate the utility of tissue-specific CRISPR screening for prioritizing candidate disease gene orthologs identified through comparative genomics platforms like Zoogle.

Bi-allelic variants of RNU4ATAC, transcribed into the minor spliceosome component U4atac snRNA, are associated to variable severity of microcephaly, growth retardation, skeletal dysplasia and immunodeficiency as main features. Previous studies highlighted the dramatic effect of U4atac deficiency on splicing of U12-type introns, which represent less than 1% of all introns in the human genome. More recently, our team evidenced a link between U4atac and the primary cilium/centrosome complex through the identification of patients carrying RNU4ATAC bi-allelic variants and exhibiting an atypical Joubert syndrome, a well-known ciliopathy. Yet, the underlying mechanisms remain elusive. Here, we further explored the link of RNU4ATAC to primary cilium and aimed at identifying ciliary U12-type intron containing genes that contribute to the brain abnormalities seen in patients. For that, we performed a transcriptomic analysis of heads of our morpholino oligonucleotide (MO)-mediated u4atac zebrafish model. Through the combined analysis of the generated dataset with those obtained from RNU4ATAC patient cells, we identified two candidate genes: TMEM107, coding for a structural protein of the cilium transition zone, and RFX7, encoding a transcription factor involved in primary cilium formation. By conducting complementary genetic approaches in zebrafish model, we showed that both gene orthologues, tmem107l and rfx7b, functionally interact with u4atac and are required for correct brain development. Altogether, our findings establish TMEM107 and RFX7 as key components of the molecular pathway linking U4atac dysfunction to ciliary defects and impaired brain development, providing new physiopathological insights and therapeutic perspectives for RNU4ATAC-related disorders.

MotivationEpigenome-wide association studies (EWAS) have identified numerous DNA methylation (DNAm) CpG sites associated with complex traits and diseases, but interpretation of those CpG sites remains challenging because in EWAS, CpGs are mostly linked to nearby genes based only on genomic proximity. Expression quantitative trait methylation (eQTM) analyses connect DNAm CpGs with statistically associated gene expression levels. However, a comprehensive, searchable resource integrating eQTMs across diverse tissues and disease contexts has been lacking. ResultsWe developed the eQTM Atlas, a web-based resource that manually curates more than 11 million DNAm-gene expression associations from eight cohorts, covering 11 tissue types, four broad disease contexts, 173,886 unique CpG probes and 20,231 unique genes. The Atlas supports gene- or CpG-searches by tissue or disease type and finding associated CpG or genes, visualization of cis- and trans-eQTMs through genome browser, heatmap interfaces across various tissues, and cohort-level data downloads. By integrating eQTM results with EWAS resources, the eQTM Atlas enables users to connect disease- or trait-associated CpGs to statistically associated genes rather than relying solely on proximity-based gene annotation, supporting functional interpretation of EWAS findings and generation of disease-specific regulatory hypotheses. Availability and implementationThe eQTM Atlas is freely available at https://shiny.crc.pitt.edu/eqtm_browser/. The web interface is implemented in R Shiny and hosted through the University of Pittsburgh Center for Research Computing (CRC). Source code is available at https://github.com/ads303/eQTM-Atlas.

Protein-protein interactions (PPIs) play a crucial role in enabling proteins to carry out their functions within various biological processes (Hui et al., 2003). Since the introduction of the yeast two-hybrid (Y2H) method for PPI detection in 1989 (Fields and Song, 1989), the identification of PPIs has become a significant focus in modern biological research. PPI goes beyond examining individual proteins, allowing researchers to establish a comprehensive network that regulates biological processes. Rice, as a key model organism in plant biological studies, has been at the forefront of PPI research. In 2008, prominent rice scientists in China called for concerted efforts to define a comprehensive protein-protein interaction network experimentally, which aimed to facilitate the prediction of the functional mechanisms operating throughout a plants lifecycle (Zhang et al., 2008). With efforts for 2 decades, the experimentally identified rice PPIs have reached over ten thousand. Several public databases have been established to systematically collate and store PPIs, including STRING (Szklarczyk et al., 2019), BioGRID (Oughtred et al., 2020), IntAct (del Toro et al., 2022), PRIN (Gu et al., 2011), RicePPINet (Liu et al., 2017) and RiceNet v2 (Lee et al., 2015). However, most PPI datasets in rice stem from computational predictions, while experiment-based rice PPI datasets are fragmented due to the lack of systematic profiling at the rice PPIome level, which largely hinders information sharing in the rice research community. To bridge this gap, we constructed the Port of Protein-Protein Interactomes (POPPIN; https://riceome.hzau.edu.cn/poppin/), an integrated database dedicated to sharing experimentally verified PPIs and functional clues in rice. Empowered by high-throughput PPIome profiling technologies and text mining assisted by a large language model (Huang et al., 2025; Liu et al., 2025), POPPIN currently has deposited over 150,451 pieces of rice PPI-related information. Additionally, POPPIN provides detailed protein information, including GO annotations, subcellular localizations, domains, trait ontology (TO) information, and hyperlinks to external biological databases. Through offering a user-friendly web interface for search and dynamic network visualization, POPPIN serves as the first large-scale, experiment-based database for searchable PPIs in rice, and has the potential to be extended to other species under this structural framework.

Limited oxygen supply or hypoxia can impair fetal development and lead to developmental disorders, but the molecular mechanism underlying this phenomenon remains poorly understood. It is also well known that hypoxia results in transcriptomic alterations and epigenetic reprogramming. Drosophila melanogaster (fruit fly) has been used for decades as a powerful model to dissect the molecular mechanisms regulating development. To better understand the role of early hypoxic stress on development, we performed single-cell joint analysis of chromatin accessibility and transcriptome to characterize the influence of hypoxia on Drosophila embryonic development. We identified hypoxia-induced alterations in both gene expression and chromatin accessibility across 22 cell groups, especially in the genes regulating organogenesis and development of neuronal, tracheal, and muscular systems, including a reduction of germ cells, suggesting a long-lasting influence of hypoxic stress at an early embryonic stage on development and reproduction. In summary, this study demonstrates that early embryonic hypoxia induces cell type- and dose-dependent changes in chromatin accessibility and gene expression, leading to distinct developmental phenotypic responses, such as reduced number of germ cells under both 3% and 5% O2. We further conclude that the tramtrack (ttk) gene is critical in germ cell development and reproduction in Drosophila melanogaster.

BackgroundMethyltransferase-like 1 (METTL1) is highly expressed in organs like the pancreas but less so in the brain. The METTL1-WDR4 complex catalyzes N7-methylguanosine (m7G) methylation in tRNA, miRNA, mRNA, and rRNA, which impacts RNA stability and function. These modifications affect mRNA translation and tRNA functionality, influencing protein production and cellular activities. Such modifications can regulate tumor growth, invasion, and metabolism by selectively controlling protein expression. MethodGene expression data from public databases were analyzed to compare METTL1 expression in normal and tumor tissues. Western blot (WB) and immunohistochemistry (IHC) were used to quantify METTL1 levels in glioma samples and assess their prognostic significance. Cell viability, migration, invasion, and proliferation were evaluated using Cell Counting Kit-8 (CCK-8), wound healing, Transwell, cell cycle analysis, and colony formation assays. RNA immunoprecipitation PCR (RIP-PCR) identified m7G methylation sites on EPHA2 mRNA, and RNA stability was assessed with actinomycin D. ResultsBioinformatics analysis revealed that METTL1 is overexpressed in gliomas, correlating with poor prognosis. Knockdown of METTL1 significantly affected cell proliferation, migration, and invasion. RNA sequencing (RNA-seq) and m7G analysis identified EPHA2 as a downstream target, influencing the cell cycle via the AKT pathway. RIP and methylated RNA immunoprecipitation (MeRIP) confirmed two m7G sites on EPHA2 mRNA regulated by METTL1. Small interfering RNA (siRNA)-mediated METTL1 knockdown in EPHA2 mutants affected mRNA stability. Rescue experiments restored cell proliferation and AKT pathway gene expression. ConclusionMETTL1 methylates EPHA2 mRNA, enhancing its stability and expression, which activates the AKT signaling pathway and influences glioma cell proliferation. METTL1 could be a potential therapeutic target in glioma treatment.

Meiosis is a fundamental process in producing both sperm and eggs, yet recombination landscapes often exhibit sexual differences, known as heterochiasmy. Since meiotic proteins are generally expressed in both sexes, the molecular mechanism driving heterochiasmy remains elusive. The -kleisin subunit gene of meiotic cohesin, Rec8, is expressed bisexually in mammals, while its putative teleost ortholog, rec8a, is expressed in a female-biased manner, presumably due to the presence of its paralog originating from the teleost-specific whole-genome duplication (TGD). Here, we elucidated the evolutionary history and expression dynamics of -kleisin genes across teleost lineages. Through comprehensive phylogenetic and synteny analyses, we revealed that major teleost lineages retain two copies of rec8 and rad21, with rec8 loci experiencing drastic chromosomal rearrangements immediately after the TGD. Using in situ hybridization and single-cell transcriptome data in medaka and zebrafish, we demonstrated a conserved sexually biased expression pattern: rec8a is predominantly female-biased, whereas rec8b exhibits male-biased expression during gametogenesis. Furthermore, comparative epigenetic analyses revealed that the conserved sexually biased expression is driven by lineage-specific cis-regulatory elements, rather than conserved ones. Motif analyses imply that regulatory rewiring by transcription factors, including foxl2l in particular, might have played a crucial role in the establishment and maintenance of this paralog divergence. Our findings highlight how whole-genome duplication and subsequent genomic and epigenetic rewiring subdivided the bisexual function of rec8, offering insights into sexually distinct meiotic regulation. HighlightsO_LITeleosts possess a unique -kleisin repertoire originating from the TGD. C_LIO_LITeleost rec8 paralogs exhibit conserved sex-biased expression during meiosis. C_LIO_LIDrastic genomic rearrangements after the duplication rewired the teleost rec8 loci. C_LIO_LIThe conserved expression pattern is governed by lineage-specific CREs. C_LIO_LIThose CREs harbor similar types of TFBSs such as Fox-family TFs. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=94 SRC="FIGDIR/small/731870v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@c8c84dorg.highwire.dtl.DTLVardef@1d65668org.highwire.dtl.DTLVardef@c2d732org.highwire.dtl.DTLVardef@1be54a2_HPS_FORMAT_FIGEXP M_FIG C_FIG

Current studies in human genomics typically rely on the standard genome reference GRCh38 which is known to be biased toward populations of European ancestry and therefore has limitations when applied to other populations. Although various graph-based pangenome references were constructed for several populations to deal with this bias, their usage in practice is currently still limited compared to linear genome references. Here we present a framework for constructing a population-specific genome reference using GRCh38 as backbone with alternate-contig awareness to enhance genomic data analysis in the target population. We demonstrated the advantages of our framework using both public and in-house Vietnamese whole-genome sequencing (WGS) datasets. Genomic variants derived from high-coverage WGS data of the 1000 Vietnamese Genomes Project (VN1K) were imported into our framework to build a Vietnamese-specific Genome Reference (VGR). VGR was then compared to GRCh38 in read alignment and variant calling using high-coverage WGS data of 99 Vietnamese individuals (KHV) from the 1000 Genomes Project (1kGP). Using Omni array genotyping data from 99 KHV samples as an independent benchmark, we found that VGR improved variant-calling precision and reduced false-positive calls compared to GRCh38. Our framework could be easily used for other populations as long as they have a variant database similar to VN1K. Our code is publicly available at github.com/VinGenome/VGR

The specificity protein 7 (SP7), also known as Osterix (OSX), is a zinc-finger transcription factor essential for osteoblast differentiation and skeletal development. Hereafter, the protein is referred to as OSX/SP7 throughout the manuscript. Although OSX/SP7 has been extensively studied in mammals and teleost fish, its developmental expression pattern in chicken (Gallus gallus) has not been described. Using an experimentally isolated chicken OSX/SP7 sequence, we examined OSX/SP7 mRNA expression during embryogenesis by in situ hybridization. As expected, OSX/SP7 expression was detected in developing skeletal elements undergoing ossification. Unexpectedly, transcripts were also observed in the neuroepithelium, retina, central nervous system, embryonic muscles and integument. Notably, OSX/SP7 expression was present in the neural tube from Hamburger and Hamilton stage 9 and persisted in the developing central nervous system until at least HH40, suggesting that OSX/SP7 functions during avian development may extend beyond osteogenesis. In addition, we reconstructed the chicken OSX/SP7 coding sequence and inferred its associated untranslated regions. The reconstructed ORF was independently supported by maternal and paternal haplotype-resolved chicken genome assemblies and retained the characteristic domain architecture of vertebrate OSX/SP7 proteins despite substantial divergence outside the DNA-binding domain. Taken together, these findings broaden the developmental landscape of OSX/SP7 expression in birds and provide new molecular resources for future studies of OSX/SP7 regulation and function during vertebrate development.

Here, we present the first graph-based Korean Pangenome Reference (K-PanRef), constructed from 14 healthy Korean individuals. K-PanRef comprises 13 high-quality diploid Korean genome assemblies (mean QV ~62.0) and KOREF1-G-TTAGGA, the first complete Korean reference genome. Integration of these assemblies generated a ~3.2-Gb pangenome graph containing ~39.3 million nodes and ~53.8 million edges, with the accumulation of common sequences (frequency [≥]10%) reaching a plateau. Additionally, K-PanRef contains ~4.3 million Korean-specific small variants and ~76.0 thousand Korean-specific SVs absent from the Chinese and human pangenome references, improving the representation of Korean genetic diversity relative to these references. To evaluate its utility for short-read-based SV analysis, we genotyped 75 whole-genome sequencing (WGS) samples, including 15 patients with early-onset myocardial infarction (MI). Although constructed entirely from healthy genomes, K-PanRef supported the identification of putative disease-relevant SVs in this exploratory application. K-PanRef-based genotyping identified ~95.6 thousand small variants and 820 SVs observed only in the early-onset MI samples. Among the early-onset MI-group SVs, 491 were absent from public databases, suggesting that they may represent previously unrecognized candidate variants related to early-onset MI. Of these, 164 SVs overlapped 134 genes, of which 89 had reported associations with 42 cardiovascular diseases or traits, including eight genes previously linked to MI. Together, these results establish K-PanRef as a valuable resource for representing Korean genetic diversity and enabling more comprehensive discovery of population-specific and novel putative disease-relevant variants from short-read sequencing data.

A protein whose removal or deletion causes significant disruption or collapse of a protein-protein interaction (PPI) network is referred to as a vulnerable protein. Such proteins may serve as valuable therapeutic or diagnostic targets in disease-associated networks. In this study, two PPI networks were constructed, one for HPV-positive and the other for HPV-negative head and neck squamous cell carcinoma (HNSCC), and the vulnerable proteins of these networks were identified by the node deletion approach. After analyzing the networks, 27 unique vulnerable proteins in HPV-positive and 72 unique vulnerable proteins in HPV-negative HNSCC were identified. Among them, one HPV-positive and seven HPV-negative HNSCC vulnerable proteins were further chosen by integrating multi-omics data. To exploit the vulnerabilities of these proteins, candidate synthetic lethal (SL) partners were predicted whose inhibition may selectively impair tumor survival. Subsequently, drug-gene interaction analysis was performed to identify inhibitors targeting the SL partners of these vulnerable proteins. Notably, in HPV-positive HNSCC, TOP2A, CHEK1, and CHEK2 genes were identified as SL partners of TTN, and their inhibitors were already clinically approved. While in HPV-negative HNSCC, ADA and MMP19 were identified as an SL partner of LMO7; TMEM45B, CDH3, and ELF3 genes were identified as an SL partner of CGN; and ZNF433 was identified as an SL partner of FLNC. However, MMP19, ZNF433, and TMEM45B inhibitors were not reported. Thus, these vulnerable proteins, including their SL partners, provide novel avenues to explore and develop more efficient and precise therapeutic and diagnostic strategies.

Despite significant advancements in microsurgical techniques in recent years, the treatment and prognosis of craniopharyngiomas remain unsatisfactory. As a central nervous system tumor located adjacent to important brain structures such as the hypothalamus-pituitary axis and accompanied by a highly inflammatory microenvironment, the tumor heterogeneity and tumor microenvironment characteristics of papillary craniopharyngiomas (PCPs) remain unclear. In this study, we integrated multimodal single-cell and spatial profiling from PCP tissue and peripheral blood mononuclear cells (PBMCs) to elucidate the tumor heterogeneity and microenvironment characteristics of PCP. Our single-cell and spatial analyses defined four specific tumor cell states in PCP, representing specific transcriptional regulatory programs and spatial heterogeneity characteristics during tumor progression. By constructing a spatial niche composed of tumor, immune, and stromal cells, we analyzed the cellular and spatial ecosystem of PCP at multiple levels to further assess the communication relationships between different tumor cell states and microenvironment cells. This study established a multidimensional molecular atlas of PCP from the perspectives of cell state, spatial structure, and microenvironment interactions, providing a foundation for understanding its biological behavior and exploring new intervention strategies.

Despite the significant technical advancement in spatial transcriptomics, its clinical usage is largely untapped. Here, we develop an integrated system, ENDO-Genome, for minimally invasive in-body transcript sampling to facilitate live spatial transcriptomic analysis of human internal organs. This is achieved by integrating a nanoarrayed biochip with existing endoscope to perform pressure-sensor-calibrated "Touch & Go" RNA extraction directly from human internal organs, including the highly vascularized liver or kidney, without the need for tissue biopsy, voiding any bleeding risks. By a demonstration using gastrointestinal endoscopy, multiplexed landscape of 55 mRNA transcripts was obtained from multiple locations of human intestinal tract via a 5-minute operation in routine examinations. Benefiting from a sequencing-free approach, each assay costs less than 10 US dollars. For the clinical study involving 15 Crohn' s disease (CD) patients, no complication case was reported out of 47 ENDO-Genome operations, showcasing the gentle deposition and excellent safety of the technique. The live spatial transcriptomics provides direct in vivo pictures of the heterogenous spatial transcriptional programs underlying CD pathological response at different intestinal locations, revealing distinct ileal phenotypes. This is manifested by unique microscale scattering of inflammation gene clusters, along with the discovery of a tissue-specific cooperative mechanisms between inflammation and RNA methylation regulations at single- or multi-cell scales.

The neural crest and neurogenic placodes, which arise from the neural plate border, give rise to morphological characteristics that distinguish vertebrates from invertebrates (Gans & Northcutt, 1983). Recent studies have suggested that embryos of ascidians, a group of tunicates that are the closest invertebrate relatives of vertebrates, possess cells that share an evolutionary origin with vertebrate neural crest cells (Abitua et al., 2012; Fatieieva et al., 2025; Ishida & Satou, 2024; Stolfi et al., 2015; Todorov et al., 2024; Waki et al., 2015) and neurogenic placode cells (Abitua et al., 2015; Ikeda et al., 2013; Liu et al., 2023; Liu & Satou, 2019; Manni et al., 2004; Mazet et al., 2005; Papadogiannis et al., 2022; Wagner & Levine, 2012). To dissect the neural plate border of ascidian embryos at the molecular level, and to gain deeper insights into evolutionary origins of the neural crest and neurogenic placodes, we comprehensively analyzed expression patterns of transcription factor genes at single-cell resolution. We demonstrated that the ascidian neural plate border consists of three domains with distinct gene expression profiles: the anterior, inner lateral, and outer lateral domains. A cross-species comparison of transcriptomes from ascidians and zebrafish suggests that the anterior domain is homologous to zebrafish neurogenic placodes, the inner lateral domain to the neural crest and tail bud and the outer lateral domain to the median fin fold ectoderm. We propose that a tripartite neural plate border was present in the last common ancestor of vertebrates and tunicates, providing a blueprint for evolutionary emergence of vertebrate morphological novelties.

RNA sequencing (RNA-seq) is widely used to investigate transcriptional programs in plant biology, yet the need to combine multiple specialized tools and bioinformatics expertise to convert raw sequencing reads into biologically interpretable results remains a major technical barrier for many plant biologists. Here, we present VizR (VIsualiZation of Rna seq), a web- based platform that integrates end-to-end RNA-seq analysis and visualization within a single integrated environment. VizR automates upstream processing, including quality control, adapter trimming, genome alignment, and transcript quantification, and connects the resulting expression data to downstream exploratory analyses. Its interface is designed to make expression patterns immediately searchable and interpretable: users can query genes through an equalizer-style expression-pattern interface, inspect expression profiles using inline heatmaps embedded in gene tables, and perform context-integrated gene ontology analysis throughout the workflow. VizR also supports comparative analysis through interactive Venn diagram module, allowing users to transfer gene sets directly from result tables. As a Docker- based application, VizR can be deployed locally and accessed through a standard web browser. By unifying automated RNA-seq processing, interactive visualization, and functional interpretation, VizR lowers the technical barrier to transcriptome analysis and provides a practical platform for plant biology research.

Cranial neural crest (CNC) cells are essential developmental contributors to the remarkable diversity of vertebrate skull shapes, yet how underlying gene regulatory networks (GRNs) evolve to produce highly derived morphologies remains a challenging question. Syngnathid fishes (seahorses, pipefishes, pipehorses, and seadragons) are an opportune family of species in which to address this problem because of their unusual and extensive cranial diversity and their loss of craniofacial patterning genes, fgf3 and fgf4. Here we investigated whether syngnathid craniofacial evolution experienced only a few localized network changes or required global rewiring of CNC GRNs. Using comparative single-cell RNA sequencing, ATAC-seq and whole genome alignments across Gulf pipefish, threespine stickleback, and zebrafish, we found that the core pharyngeal arch CNC gene network is notably conserved in syngnathids despite their derived morphology. However, we identified key local changes including expression of fgf22 in percomorph CNC-derived pharyngeal arch cells that is not shared with more basally diverging zebrafish, as well as syngnathid-specific changes in conserved regulatory elements associated with the genes ece1 and spry2. We propose that, while loss of fgf3 expression causes severe craniofacial defects in zebrafish, pharyngeal CNC expression of fgf22 in the percomorph fish lineage provided functional redundancy and relaxed constraint on fgf3/4, and that altered regulation of Fgf pathway modulators could contribute to craniofacial elaboration. Our findings support a model in which local GRN modifications, rather than widespread network rewiring, underlie the evolution of derived syngnathid craniofacial structures. SignificanceUnderstanding how developmental genetic changes drive the evolution of unique traits remains a long-standing challenge in biology. In the case of syngnathid fishes (pipefishes, seahorses, and seadragons), previous genomic studies identified candidate craniofacial gene losses which are proposed to relate to their elongate and derived heads, but the developmental impact of these losses is unknown. Through gene expression and comparative genomics analyses, we find that these fishes have distinct changes to craniofacial gene regulatory networks including gene content losses and gene expression gains and losses. Our study suggests that morphological adaptations may arise from multiple key changes within largely conserved developmental regulatory networks.

Oaks have been called the "worst-case scenario" for the biological species concept because of the presumed importance of introgression among species. However, we still lack an estimation of the importance of genome-wide introgression across the eight sections of the Quercus phylogeny. Nor do we have a complete characterization of what is introgressed and how it functionally influences the recipient species. Firstly, we used 168 whole genomes from 67 Quercus species belonging to the eight sections of the Quercus phylogeny and a super pan-genome based on 9 species to generate a phylogeny and test for incomplete lineage sorting and hybridization. Introgression was widespread and occurred within and across the eight sections. Ancient introgression occurred between sections, while more recent introgression took place within sections. Secondly, using genome-wide transcriptome and methylome from a subset of 19 species sampled from four sections we demonstrated that introgressed genes tend to be more highly and uniformly expressed than non-introgressed ones. This difference in expression level between introgressed and non-introgressed genes mostly results from cis-regulation, and both methylation and chromatin accessibility explain the high and conserved expression of introgressed genes across species. In summary, introgression has been a major evolutionary force during the whole evolutionary history of Quercus and had a significant impact on fundamental processes such as gene regulation and expression. Significance StatementMany oak species have been known to exchange genes during their evolution. However, there are more than 400 oak species today, divided into eight phylogenetic sections and distributed over four continents, and the extent to which introgression has affected the whole oak phylogeny remains unknown. Using whole genomes from 67 oak species belonging to the eight sections we show that introgression occurred within and across the eight sections. Ancient introgression occurred between sections, but more recent introgression took place within section. Using genome-wide gene expression and methylation data from a subset of species we also demonstrated that introgression affected gene regulation and expression: introgressed genes are more highly and uniformly expressed than non-introgressed ones.

A central task in population genetics is to identify genetic diversity in a population containing a number of individuals. In recent years, with the development of the third generation sequencing (TGS) technology, pan-genome research has become a hot topic. Although graphical representation has been a popular way to represent the pangenome, few works have attempted to describe it in a more mathematical way. In this paper, we used 79 high-quality assembly data of third-generation sequencing in yeast (includingSaccharomyces cerevisiae and Saccharomyces paradoxus) to construct the graph pangenome, and introduced the Primitive GLMY (Grigoryan-Lin-Muranov-Yau) Homology in algebraic topology to quantitatively represent the pan-genome. We further made an intriguing attempt to conduct a population genetic analysis of this resulting dataset from the topological features of the graph pangenome. We found that there was good agreement between the obtained results and the biological context. We believe this study has developed a method for population genetic analysis of the genetic diversity of genome structural variation.